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Abstract 

We present a search for the decays —>■ e^e~ , — >■ jji^ , and B^ e^/U^ in data collected at 
the T{AS) with the BABAR detector at the SLAC B Factory. Using a data set of 54.4 fb"\ we find 
no evidence for a signal and set the following preliminary upper limits at the 90% confidence level: 
B{B^ e+e-) < 3.3 x lO^^ B[B^ ^i+fi-) < 2.0 x 10"^ and B{B^ e=t/x=F) < 2.1 x lO"'^. 
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1 Introduction 




Figure 1: Feynman diagrams for l+l- in the Standard Model. 

since they involve a 6 — > d transition and require an internal quark annihilation within the B meson 
which further suppresses the decay by a factor of {Jb/Mb)'^ 2 x 10~^ relative to the electroweak 
"penguin" b dj decay. In addition, the decays are helicity suppressed by factors of {mi/mB)'^ ■ 
B^ decays to leptons of two different flavors, B^ — s- e^/z^, violate lepton flavor conservation and are 
therefore strictly forbidden in the SM, although permitted in extensions to the SM with non-zero 
neutrino mass 0]. The B^ e^fi'f channel includes both B° e+//- and B^ 6"//+. To date, 
B'^ — > £~^£~ decays have not been observed. The current best limits from the CLEO ||2| and Belle ||3| 
collaborations are compared with the SM expectations in Table |l[ 

B^ decays involving tau leptons involve either a soft electron, muon, or meson and 

missing energy (from one or more neutrinos), and require a rather different search strategy than 
that presented here for the B^ — > e~^e~, B^ and B^ — > e^/x^ channels. The presence 

of the soft electron or muon in the tau channels eliminates them as a source of background to the 
non-tau channels. 

Since these processes are highly suppressed in the SM, they are potentially sensitive probes 
of physics beyond the SM. Although the branching ratios for these processes are not significantly 
enhanced in the Minimal Supersymmetric Standard Model (MSSM), various non-minimal super- 
symmetric models predict branching ratios that are significantly larger than those of the SM. Also, 
multi-Higgs-doublet models with natural flavor conservation have extra charged Higgs particles 
which replace the SM l^-boson in the box diagram of Fig. [l| and may enhance the B^ — > 
branching ratios by an order of magnitude . Similarly, in models with an extra vector-like down- 
type quark [^, flavor changing neutral currents (FCNC) involving the boson are induced due to 
the different isospin charge of the exotic quark. These models predict the rate for B^ — s- £~^i~ to be 
about two orders of magnitude larger than the expected SM rate In addition, B^ — > ("^t' de- 
cays are allowed in specific models containing leptoquarks Q and supersymmetric (SUSY) models 
without i?-parity conservation. Furthermore, flavor violating channels such as B^ — > e^/x^ could 
be enhanced by leptoquarks or ii-parity violating operators in SUSY models. 
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As shown in Table sensitivity even to models which produce an enhancement of a few orders 
of magnitude to the SM rates for these rare decays is beyond current experimental capabilities. 
Observation of a i^£~ decay would, in consequence, provide clear evidence for physics beyond 

the Standard Model. 

Table 1: The expected branching ratios in the Standard Model and the current best upper 
limits (U.L.) in units of lO"'^ at the 90% C.L. from CLEO [|] and Belle g. In addition, the 
measured number of events A'obs in the signal region, the expected background in the signal region 
A^^p, and the reconstruction efficiency e are quoted. CLEO's analysis was performed on a data set 
corresponding to 9.1 fb~^. Belle's measurement was performed on a data set of 21.3 fb~^. 



Decay 
Mode 


SM Expect. 


CLEO 


Belle 


U.L. 


Nohs 


cxp 


s[%] 


U.L. 


A^obs 


cxp 


s[%] 




1.9 X 10-^^ 


8.3 





0.11 ±0.07 


31.1 ±0.4 ±2.4 


6.3 


1 


0.6 ±0.8 


31.3 ±2.4 




8.0 X 10-" 


6.1 





0.22 ±0.07 


42.4 ±0.5 ±3.2 


2.8 





0.7 ±0.8 


40.0 ±4.3 






15 


2 


0.49 ±0.20 


43.6 ±0.5 ±7.1 


9.4 


3 


0.7 ±0.8 


35.8 ±3.2 



2 The BaBAR detector and dataset 

The data used in these analyses were collected with the BABAR detector at the PEP-II e~^e~ 
storage ring during the years 2000 and 2001. The sample corresponds to an integrated luminosity 
of 54.4 fb~^ accumulated on the T{4S) resonance ( "on-resonance" ) and 6.4 fb~^ accumulated at 
a center-of-mass (CM) energy about 40MeV below the T{4S) resonance ("off-resonance"), which 
are used for the non-resonant qq background studies. The on-resonance sample corresponds to 
(59.9 ± 0.7) X 10^ BB pairs. The collider is operated with asymmetric beam energies, producing a 
boost (/37 = 0.55) of the T(45') along the collision axis. 

BABAR is a solenoidal detector optimized for the asymmetric beam configuration at PEP-II and 
is described in detail in Ref. |^]. The 1.5 T superconducting magnet, whose cylindrical volume 
is 1.4 m in radius and 3 m long, contains a charged-particle tracking system, a Cherenkov detec- 
tor (DIRC) dedicated to charged particle identification, and a central electromagnetic calorimeter 
consisting of 5760 CsI(Tl) crystals. A forward endcap electromagnetic calorimeter consists of 820 
CsI(Tl) crystals. The segmented flux return, including endcaps, is instrumented with resistive plate 
chambers for muon and identification. This subsystem is referred to as the instrumented flux 
return (IFR). The tracking system consists of a 5-layer double-sided silicon vertex tracker and a 
40-layer drift chamber filled with a gas mixture of helium and isobutane. 

3 Analysis method 

The presence of two charged high-momentum leptons provides for a very clean signature for the 
three decay modes under consideration. In the CM we require two oppositely-charged high- 
momentum leptons (i.e. \p£\ ~ 771^/2) from a common vertex consistent with the decay of a B^ 
mesonQ. Since the signal events contain two B^ mesons and no additional particles, the energy of 

^Charge-conjugation is implied throughout this paper. 
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each in the CM frame must be equal to the or e beam energy. We therefore define 

"^Es = ^{EL.j'-iT.pty (1) 

AE = E\A? + ^-Keam, (2) 
i 

where -E-beam beam energy in the T(4S) CM frame, p* is the CM momentum of particle i 

in the candidate S'^-meson system, and rrii is the mass of particle i. For signal events, the beam- 
energy-substituted mass, m-ESj peaks at m^. The quantity AE is used to determine whether a 
candidate system of particles has total energy consistent with the beam energy in the CM frame. 
We require that the beam-energy substituted mass, tues^ be very close to the mass of the B^ meson 
and that AE be close to zero 

To remove background from lepton misidentification, we require tight electron and muon iden- 
tification criteria. The electron identification relies on E/p, the ratio of energy deposited in the 
calorimeter to the momentum of the particle at the origin, the lateral and azimuthal shower profiles, 
and the consistency of DIRC Cherenkov angle with an electron hypothesis. The muon identifica- 
tion relies on the total number of hits in the IFR, the distribution of hits in the different layers, 
the amount of energy released in the calorimeter, and the number of interaction lengths which the 
track has traversed. Suppression of background from non-resonant qq production is provided by a 
series of topological requirements. In particular, we employ restrictions on the overall magnitude 
of the event thrust and on the magnitude of the cosine of the thrust angle, 6t, defined as the angle 
between the thrust axis of the particles that form the reconstructed B^ candidate and the thrust 
axis of the remaining tracks and neutral clusters in the event. We also cut on the total multiplicity 
of both charged tracks and neutral particles by means of the variable A^muit defined as 

iVmult = iVtrk + ^N^, (3) 

where A'trk is the total number of tracks in the event and is the number of photons found with 
an energy E^ > 80 MeV. We require A^^muit ^ 6.0. This variable is especially useful in the rejection 
of radiative Bhabha events. We also require that the total energy in the event be less than 11 GeV. 
This cut is effective in reducing background from two photon events. 

B^ i+i- candidates are selected by simultaneous requirements on the energy difference 
AE and the energy-substituted mass niES- The size of this "signal box" is chosen to be roughly 
[-1-2, — 2](T of the expected resolution in AE and [-1-2, —2]a for niEs, optimized for the best upper 
limit. In the cases of the B^ — > e^e" and B^ —>■ e^/j,^ decay modes, the signal box sizes in 
AE are relaxed to roughly [+2, —3]a and [-1-2, —2.5]a, respectively, to accommodate a tail in the 
distribution resulting from final state radiation and bremsstrahlung. Table ^ gives the niEs and 
AE resolutions for the three signal channels. Figure § illustrates the tues and AE distributions 
for Monte Carlo (MC) signal events. For the AE distribution, the tail due to final state radiation 
and bremsstrahlung is well described by a "Novosibirsk" function jlOI J. Table ^ summarizes the 
ruES and AE cut values used to define the different boxes in the {AE,mEs) plane. 

We also chose to optimize the cuts on the magnitude of the overall event thrust |r[ and | cos 9t\ 
simultaneously (due to the large correlation). The optimal selection criteria for all three channels 
was found to be [cos 0^1 < 0.84 and |T[ < 0.9. Figure ^ illustrates the distributions of the 
multiplicity and event shape variables in signal and background MC, which are dominated by 
non-resonant cc- and nds-continuum processes, but also include small components from bb and 
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Figure 2: Reconstruction of B meson candidates in e^e~ MC (top row), B^ /U"*"^" MC 

(middle row), and B^ — > e^/i^ MC (bottom row) with the beam-energy substituted mass uies 
(a) and A£^, the difference between the beam energy and the energy of the i3-meson candidate 
in the CM frame (b). A "Novosibirsk" function [0] is used to obtain the widths of the core 
distribution. Figures (c) show the distribution of l\E vs tues- The smaher box on the right defines 
the signal box. The tail visible on the lower side of the signal box is due to final state radiation 
and bremsstrahlung. The larger box defines the Grand Sideband (GSB) region. 
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Table 2: Definition of the signal and Grand Sideband boxes in the {AE,mEs) plane. The units for 
niES are [GeV/c^] and the units for AE are [GeV]. 









B" e±M=F 


Box Name 




AE 




AE 


rriES 


AE 


Signal Box 
Grand Sideband 


5.273-5.285 
5.200-5.260 


-0.105-0.050 
-0.400-0.400 


5.274-5.285 
5.200-5.260 


-0.050-0.050 
-0.400-0.400 


5.274-5.284 
5.200-5.260 


-0.070-0.050 
-0.400-0.400 



r events. All selection criteria have been applied except for the cut on the variable illustrated. 
The efficiencies of the full selection are given in Table The systematic error on the efficiency is 
determined by a comparison of the control sample — s- J/tjjKg, with J/ijj e^e~ for B^ — > e'^e~ 
and J/il^ ^ IJ-~^f^^ for B^ respectively in data and MC simulation. These comparisions 

found the dominant uncertainty on the signal efficiency to be the resolution and scale of AE, 
contributing 4.4% and 2.6% for the B^ — > e~^e~ and B^ channels respectively. Since there 

is no appropriate control sample for the B^ — > e^/i+ channel, we conservatively set it to be equal to 
the largest error obtained from the systematic study for the other channels. This yields systematic 
errors for all of the main cuts except the multiplicity cut. The systematic error associated with the 
remaining cuts is determined with a comparison of signal MC samples based on GEANT3 |11] and 
GEANT4 |]l^]. These two simulations employ different material models where the latter is considered 
to be more accurate. 
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Figure 3: Distributions of the main cuts in signal MC and background MC. The histograms are nor- 
malized to equal area. The cut values are indicated by arrows. The composition of the background 
MC is explained in the text. 



Table 3: niEs and AE resolutions for the three signal channels. 





B^ e+e- 


B" fi+fi- 




a{mEs) [MeV/c^] 
a{AE) [MeV] 


3.0 ±0.2 
29.3 ±0.9 


2.6 ±0.1 
24.7 ±0.3 


2.7 ±0.1 
26.8 ± 0.4 



We estimate the background level in the signal box from the data. The background expectation 
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is dominated by different sources in the tliree cliannels. For tfie — > e^e^ cliannel, the background 
expectation is dominated by pairs of true electrons from cc events and two-photon events. For the 
fi^ fi^ channel, about 50% of the total background is due to misidentified hadrons (in 
combination with a real muon). Two-photon processes do not contribute to the background for 
this channel. For the B^ — > e^/i^ channel, the background is composed of real electrons and fake 
muons. Two-photon processes contribute strongly to the background. 

To compare background distributions in data and Monte Carlo event samples, we use the "Grand 
Sideband" (GSB) box as defined in Table This box is also used to estimate the functional behavior 
of the /S.E dependence of the background. 

We estimate the background in the signal box assuming that it is described by the ARGUS 



function |13| in m^;^ and an exponential function in Ai?. The shape parameters of the functions are 
determined in two different ways: from fitting the data sidebands and from fitting a high statistics 
fast Monte Carlo cc sample. In both cases, we determine the normalization in the grand sideband 
and derive the number, A^'sigBox, of expected background as follows: 

/gsb 



AT J.igBox-'^'"-^^^"-"-'^'^ ^ JsigBox ^^"""^ "^"""^ ^ Ar 



where Nq^^ is the number of background events found in the GSB, and / and g are the shapes as 
determined by the ARGUS and exponential fits. The total background expectations for the three 
channels are given in Table ^. 

The actual contents of the signal box were not revealed until the selection criteria and systematic 
error estimates were frozen. This technique, often referred to as a blind analysis, is adopted to avoid 
possible experimenter bias. 



4 Results 

When the contents of the signal box were revealed, one event was found in the B^ e'^e~ 
channel and no events were found in the other channels as summarized in Table ^ The {rriEs, AS) 
distributions from data for the three channels are shown in Fig. ^ As can be seen from Table ^, 
the number of events found in the signal box are compatible with the expected background. 

We do not perform background subtraction for the determination of the upper limit of the 
branching fraction. Not yet accounting for the systematic uncertainties, the upper limit on the 
branching fraction ;BuL(stat) is calculated as ;BuL(stat) = Nul/S = Nul/{£ x {N^o + N^o)), 
where Njjl is the upper limit on the number of observed events Aobs and S is the sensitivity. 
N^o + N^o is equal to the number of T(4S') decays, since we assume equal production of B'^ and 
B^ in r{4S) decays. 

In order to include our systematic uncertainty in the determination of the upper limit, we follow 



the prescription given in Ref. |14]. Assuming a normal distribution for the uncertainty in 1/5", the 
systematic uncertainty is accounted for by convolving the Poisson probability distribution for the 
assumed branching fraction with a Gaussian error distribution for 1/5". The systematic uncertainty 
on the signal efficiency is found to be 8.2% for the B^ e^e~ channel, where the main contribution 
comes from the modeling of the mEs and AE resolutions, and from the uncertainty in the efficiency 
of the multiplicity A^muit requirement. For the B^ — > iX^ channel, the systematic uncertainty 
on the signal efficiency was found to be 4.7%, with the primary contribution again being from the 
modeling of the niEs and AE resolutions. For the B^ — > e^fi^ channel, the systematic uncertainty 
on the signal efficiency was taken to be the same as for the B^ — > e'''e~ channel. The total 
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Figure 4: Distributions from data of niEs vs for 
right), and e^jX^ (bottom). 



e+e- (top left), B^ 
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Table 4: Summary of the analyses. A'exp is the number of expected signal events, assuming a 
branching fraction of 1.9 x 10~^^ for the — > e"*~e~ channel and 8.0 x 10~^^ for the — > ^jl^ ^~ 
channel. A'obs is the number of observed events in the signal box. is the expected number of 

background events in the signal box. B\]\^{B^ t^l~) is the upper limit on the branching ratio at 
the 90% C.L. 



channel 


iVcxp 


iVobs 


cxp 




Sul(B" ^^+^') 


e+e- 


1.9 X 10-« 


1 


0.60 ±0.24 


19.3±0.4,tat±1.6sys 


3.3 X 10-^ 


B° ^ 11+ ii- 


3.2 X 10-2 





0.49 ±0.19 


18.8 ± O.Sstat ± 2.0sys 


2.0 X lO"'^ 


^ e+^JL+ 







0.51 ±0.17 


18.3±0.4stat±1.5sys 


2.1 X 10""^ 



uncertainty was calculated by summing in quadrature the systematic uncertainty on the number 
of BB events, the signal efficiency, and the statistical error on the signal efficiency. 

As summarized in Table ^, the resulting preliminary upper limits for B{B^ e+e"), B{B^ 
fi-^^jL'), and B{B^ e^f^^) are 3.3 x 10"^, 2.0 x 10"^, and 2.1 x 10"^ respectively. 
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